Discrete versus continuous wires on quantum networks 
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Mesoscopic systems and large molecules are often modeled by graphs of one-dimensional wires, 
connected at vertices. In this paper we discuss the solutions of the Schrodinger equation on such 
graphs, which have been named "quantum networks". Such solutions are needed for finding the 
energy spectrum of single electrons on such finite systems or for finding the transmission of electrons 
between leads which connect such systems to reservoirs. Specifically, we compare two common 
approaches. In the "continuum" approach, one solves the one- dimensional Schrodinger equation on 
each continuous wire, and then uses the Neumann-Kirchoff-de Gennes matching conditions at the 
vertices. Alternatively, one replaces each wire by a finite number of "atoms", and then uses the 
tight binding model for the solution. Here we show that these approaches cannot generally give the 
same results, except for special energies. Even in the limit of vanishing lattice constant, the two 
approaches coincide only if the tight binding parameters obey very special relations. The different 
consequences of the two approaches are demonstrated via the example of a T-shaped scatterer. 



I. Introduction 

The quantum mechanics of electrons on networks, 
made of one-dimensional (ID) wires which are con- 
nected at vertices, has been the subject of active re- 
search for a long time. This research probably started 
with the study of organic molecules, in 1953^!^ How- 
ever, it gained a large momentum in 1981, when P. G. 
de Gennes considered the linearized Landau-Ginzburg 
equations on networks, in connection with percolating 
superconductorsi^ This work was later extended by 
Alexanderi^ In the last 25 years, quantum networks 
have been used as theoretical models of mesoscopic sys- 
tems, which have become experimentally available due to 
the advancements in microfabrication. In these models, 
the ID wires represent limits of three dimensional wave 
guides, which are assumed to be sufficiently narrow to 
justify keeping only the lowest transverse mode. Quan- 
tum networks are also used to model novel molecular 
devices, where each wire consists of discrete atoms. 

In this paper we compare several approaches for solv- 
ing the single electron Schrodinger equation on such net- 
works. Basically, this solution consists of two stages. 
First, one solves the equation on each wire. Denoting 
the vertices at the ends of the wire by a and /3, one can 
express the solution inside the wire, ipafiixap) (where 
Xap = at vertex a and Xap = La/3 at vertex /3) in 
terms of the wave functions at its ends, ^^ and ^^3. At 
the next stage, one must employ some matching condi- 
tions at each vertex. Using these matching conditions, 
and eliminating the wave functions inside each wire, one 
ends up with a set of discrete linear equations for the 
wave functions at the vertices. 
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where the coefficients Map are determined by the ID 
solution on each wire and by the matching conditions. 
Obviously, the detailed solution for the network depends 



on the latter. 

There have been several approaches to the issue of the 
matching conditions. Some of the earlier papers in meso- 
scopics used the scattering approach i^i^ For a vertex 
with Ma wires, this approach constructs an Ma x Ma 
unitary scattering matrix;^ which relates the amplitudes 
of the Ma incoming (plane) waves to the Ma outgoing 
waves. Given this matrix, one can solve for the ampli- 
tudes of the two waves in each wire, and thus find a full 
solution for the wave functions on the network. For the 
examples treated in these early papers, the scattering 
matrix was parametrized arbitrarily. To find explicit val- 
ues for the elements of the scattering matrix, one must 
solve a more "microscopic" model, as done e.g. in Refs. 
113 and nil. 

A widely used second approach follows dc Gennes ^1^ 
and uses the so called Kirchoff (or Neumann) matching 
condition at the vertices . ^^'"'^'^i-'^^'-'^^i-'^^'-'^'^i"'^^ Given the con- 
tinuous solutions Tpap (xap) OH each of the Ma wires which 
meet at vertex a, one requires continuity of the wave 
function, -00,^(0) = ^q, and also the condition 
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where the sum is over the Ma wires which meet at vertex 
a. Since the outgoing probability current on wire (a/3) 
is proportional to 



lafi = 3[(7/'a/3)*(a-0a/3/l9Xa/3)], 



(3) 



Eq. ^ supplies a sufficient condition for the conserva- 
tion of the probability at this vertex. 

However, Eq. ([2|) docs not represent a necessary con- 
dition for the conservation of current. A more general 
condition is given by 
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with any real Aq, ,i9i20i21,22 Multiplying both sides of this 
equation by (^a)* yields a real number, which does not 
contribute to the current in Eq. ([3]). For a vertex with 
two wires, Aq represents the strength of a delta- function 
potential at the vertexi^ An infinite value of Aq in Eq. 
([4]) implies that '^a = 0, which is known as the Dirichlet 
matching condition. In this limit, there is no probability 
current through this vertex. Thus, finite values of A^ in- 
terpolate between the Neumann and the Dirichlet limits. 
However, for real networks one still needs to derive A^ 
from some microscopic equations. 

An alternative approach, which avoids some of these 
ambiguities, is based on the tight binding model, as done 
e.g. in Ref. [l3- In this approach, the system is described 
by a set of localized wave functions, associated with dis- 
crete "atomic" sites {n}. This approach is often used for 
molecules, where it is equivalent to the method of linear 
combination of atomic orbitals (LCAO). It is also used 
for the calculation of energy bands in periodic lattices;^ 
The discrete Schrodinger equations for the amplitudes 
{ipn} of the "atomic" stationary wave functions for an 
energy e are then given by 



(e - en)ipn 



/ ^ ^nmV^n 



(5) 



Here, e^ and —Jnm are the diagonal and the off-diagonal 
matrix elements of the Hamiltonian, in the basis of the 
atomic wave functions. For real molecules or crystals, 
these matrix elements can be calculated from the atomic 
wave functions and from the microscopic potentials seen 
by the electron i^ The tight binding equations have also 
been applied to a variety of specific mesoscopic quantum 
networks, where the "atoms" represent quantum dots; 
see e. g. Refs. [13, EilHIlll^ and references therein. 

For ID chains, the results of the discrete tight-binding 
model approach those of the continuous Schrodinger 
equation when one sends the lattice constant of the for- 
mer to zero (see below). However, in this paper we show 
that generally this limit is problematic for more complex 
networks. Section II presents the derivation of Eq. ([T|) 
for continuous wires and for discrete wires (described by 
the tight binding Hamiltonian), and contains a critical 
comparison between them. Section III demonstrates the 
difference between the two approaches for the simplest 
example of a T-shaped scatterer, and Sec. IV contains 
our conclusions. 



II. The equations for the vertex wave functions 

As stated above, applying the Neumann matching con- 
ditions to the wave functions describing the continuous 
wires yields Eq. ^. Here we first rederive this equation 
for the more general matching condition (|3|). then de- 
rive a similar equation from the solution of the discrete 
model, and finally compare the two resulting equations. 



A. Continuous graphs 

For simplicity, we follow much of the literature and as- 
sume that there is no potential energy within the wires. 
We start by solving the continuous ID Schrodinger equa- 
tion for a free particle on each ID wire. For a given en- 
ergy, e, the solutions are linear combinations of e*'^^ and 
g-jfcx^ with £ = h^k'^ /{2m). Here we used the simplify- 
ing assumption that the effective mass m is the same on 
all the wires (in general one could have a different wave 
vector k for each wire, for the same energy e). We next 
express the solution ijjap on the wire {af3) in terms of the 
wave functions at its two ends j^^i^^ which we denote by 
^a and */3: 



Ipa/liXaff) = *« 
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Substituting these expressions into the matching condi- 
tion (|4]) then yields 



^fc[-*QCot(/cLQ^)-t-*^/sin(fcLQ^) 

/3 



(7) 



Generally k could have different values for different wires. 
Assuming the same k for all wires, this can be rewritten 
in the form of Eq. ([1]), with 
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Indeed; Aq — > oo implies that A/qq —^ cx), and therefore 
^Q = 0, as in the Dirichlet limit. The solution of the set 
of linear equations ^ then either yields ratios among 
the ^'s at a given energy e (for a scattering problem) or 
gives the set of spectrum of allowed e's (for an eigenvalue 
problem). 

B. The tight binding model 

We now replace the {a(3) wire by a discrete chain 
of "atoms". For simplicity, we assume that all the 
wires have the same lattice constant a. The wire {a(3) 
now has Nap + 1 "atoms" , at sites Xap = na, where 
n = 0, 1, 2, ..., Nap- The length of the wire is then given 
by Lap = aNafj- Assuming no potential energy within 
each wire also allows us to set the diagonal site energies 
e„ within each wire equal to zero. For simplicity we also 
assume uniform nearest-neighbor "hopping" matrix ele- 
ments, Jn,n±i = J- For < n < Nap, the tight binding 
equation ([5]) becomes 

f^ipapin) = -Jiipapin- 1) +i^ap{n+l)], (9) 

with solutions e'*^"" or e"**^"" and with the dispersion 
relation 



e = — 2Jcos(fca). 



(10) 



(Again, if different wires have different J's or different 
a's then each wire will have its own k, for the same e). 



Fixing the wave functions at the ends of the wire, 
V'q/3(0) = ^Q and ^apiNajj) = ^/3, the solution along 
the chain (a/?) becomes 

sin[fca(JVa/3 - n)] sin(fcQn) 

tpaf3[n) ^ . . ^ ^a + . ,, ,. r*/3- (H) 



sin(fca-/VQ,^) 



sin(A:a7VQ,^) 



As might be expected, this solution is just a discrete ver- 
sion of Eq. (0, but with the modified dispersion relation 

We next write down the tight binding equation ([5]) for 
the vertex site a, where we do allow for a non-zero site 
energy ea- 



(e- ea)'^a = -Jy-ipapin = 1). 



,3 



(12) 



Substituting Eq. ^^, and using Eq. UHl) (with the 
same k on all wires), this again yields Eq. ([1]), and Mafj is 
again given by Eq. ([8]). However, the diagonal coefficient 
is different: 

Maa = yj[cot(fcaA^Q^) — cot(fca)] 

+ 2cot(fca) -f-eQ/[Jsin(fca)]. (13) 

Note that in some sense Eq. ([T]) has the same form as 
the discrete Schrodinger equation (O for the vertices of 
the network. However, unlike the Jnm's and e„'s of that 
equation, the coefficients Maa and Maf} do depend on 
the energy e which characterizes the solution. Note also 
that these renormalized coefficients reduce back to those 
in Eq. ([5]) when all the wires reduce to single bonds. 

Nap = 1. 

C. Comparison 

Comparing Eq. (|13p with Eq. ([5]) , we see that the two 
equations would become identical only when 

Xa/k = ea/[Jsm{ka)] - {Ma - 2)cot(A:a). (14) 

(Again: the terms which depend on ka should be replaced 
by sums of such terms with different fc's if the wires have 
different parameters). The on-site energy on the vertex 
Eq represents a diagonal matrix element of the Hamilto- 
nian, and therefore it should not depend on the energy e 
which represents the specific solution for the Schrodinger 
equation. Therefore, Eq. (J14p implies that the "equiva- 
lent" continuum Hamiltonian would have parameters Aq 
which generally depend on e. However, the A^'s represent 
the matching conditions, and these are also characteris- 
tic of the system and not of its solutions. Therefore, 
in general one cannot choose energy-independent vertex 
parameters for which the two approaches yield the same 
solutions everywhere. 

Interestingly, the difference between the two ap- 
proaches is largest when cot(fca) is large, namely near the 
band edges. Usually, divergences near the ID band edges 



are associated with van Hove singularities, which cause 
the divergence of the density of states. For this reason, 
the use of tight binding equations is often restricted to 
the vicinity of the band center, i. e. ka ~ tt/2. However, 
even in this region the two models become identical only 
at the band center 6 = and not at any other energy. 
Furthermore, although one can choose the electron en- 
ergy (i. e. the Fermi energy) in a scattering calculation, 
where the spectrum is continuous, this choice may be 
limited in a calculation of the whole (discrete) spectrum 
of a finite network, or in a calculation which may involve 
several bands (or discrete energies) with gaps between 
them. 

It is interesting to consider this comparison for a ver- 
tex on a linear chain, when Afa — 2. Even in this case, 
the parameters Aq or e^ would have to depend on energy, 
except for the trivial case Aq = Cq, = 0. The latter condi- 
tion simply implies that this vertex is equivalent to any 
other point along the wire on which it "sits" . Although 
appropriate in some cases, this cannot be generally valid. 
For example, this does not suffice to describe the config- 
uration of two wires which meet at a finite angleJ^ Such 
a break in the straight line should cause scattering even 
if the wire itself is free of any scattering centers. To de- 
scribe this scattering, one would need to introduce non- 
zero values for £„ and/or Aq. At this stage, we are not 
aware of a systematic study of this issue. 

A priori, one might have thought that the two ap- 
proaches should coincide when the lattice constant in the 
tight binding model approaches zero. Indeed, as noted 
above, the two models yield the same answers within each 
(uniform) linear chain. In the limit of vanishing a, Eq. 
^ can be written as 



{e + 2J)xJja(3in) 



- -J[lpal3in - 1) + i>al3{n + 1) 
2Va/3(«)J ^ ~Ja 



dx^ 



(15) 



This reproduces the continuum Schrodinger equation, 
provided we move the zero of energy to the bottom of 
the band, — 2J, and identify the effective mass m via 
fi^ /(2m) ~ Ja^ , requiring that either tti ^ cxd or J ^ oo 
and a ^ 0, such that Ja^ is kept ffirite. 

We next turn to the Schrodinger equation at the ver- 
tex, Eq. UH). In the limit of small a, we rewrite this 
equation as 



(e-eQ+7VQj)*Q 






= -J> jV'a/3(l)-^./3(0)] 



-^«E 



dlpal, 



dx 



Q/3 



(16) 



For a finite effective mass m, (e -|- 2 J) ^ Ja^k^ remains 
finite when a ^ 0, and this equation becomes equivalent 
to Eq. dH) with 



Ao 



(A/; - 2) J]/( Ja). 



(17) 



In fact, this equation also follows from Eq. P^ in the 
limit of vanishing a. For Afa 7^ 2. the two models coincide 
in this limit only if either A„ ^ 00 (i.e. the Dirichlct 
matching condition which implies ^a = 0) or Ca ~ (A/^a — 
2) J - l/a2 



00 1. e. 



(A^a-2)J + M„/a, 



0. 



(18) 



The correction term, /Xq,, determines Aq. via Aq = 
Ha /{'Jo?')- To reproduce the Kirchoff-Neumann-de 
Gennes matching conditions ((2), one would need to have 
^la = 0. 

At this stage, we have no direct information on how the 
tight binding Hamiltonian matrix elements vary with a 
for small a. In fact, the tight binding model was con- 
structed to deal with small overlaps between the local- 
ized wave functions, which is no longer valid when a ^ 0. 
However, it would still be interesting to have explicit mi- 
croscopic mathematical calculations of the tight binding 
matrix elements, and to check whether indeed they all di- 
verge as 1/a?. It would be even more interesting to find 
out whether the limiting tight binding matrix elements 
obey the relation P^ . 

An exception arises for a vertex on a linear chain. Ma = 
2. In this case, Eq. pT)) reduces to Aq '^ ea/{Ja). For 
a finite effective mass on the wires, Ja^ remains finite 
when a — > 0, and therefore Aq — > whenever e^ is finite, 
reproducing the Neumann matching condition. 

Some of the papers which use tight binding models 
to describe mesoscopic systems actually use e„ = for 
all the sites, except for a few resonant sites for which 
e„ is varied by an external gate voltage. As explained, 
such a variation does not affect the Neumann matching 
conditions for A/'q = 2 when a —^ 0. However, if such 
gate voltages are applied to a vertex with more than two 
wires, then the parameter Aq for that vertex will increase 
with decreasing a, crossing over to the Dirichlet condition 
and to the vanishing of the wave function on that vertex 
when a — > 0. 



III. The T-shape scatterer 

Quantum networks have been studied for a variety of 
geometries. Here we demonstrate the differences between 
the discrete tight binding and the continuum models on 
one of the simplest geometries. We calculate the trans- 
mission of an electron passing through a T-shaped device, 
consisting of a main wire attached to a stub of length L 
perpendicular to the wire at the origin x = (Fig. [1]) J^ 
Assuming that the electrons come from the left hand side, 
we expect a reflected wave on the left and a transmitted 
wave on the right hand side. 

In the continuum approach, the wave function on the 
main wire is 



^«,(2^) = e''=^+re-''=^^ x < 0, 



^^ix)=te''^ 



a; > 0, 



(19) 



FIG. 1: A wire with a stub. 

where r and t are the reflection and transmission ampli- 
tudes. Thus, ipui{0) = 1 + r = t. From Eq. ^, the wave 
function on the stub is given by 



ips{y) = ipwiO) 



sm[k{L - y)] 
sm{kL) 



(20) 



vanishing at y = L. Using Eq. ^ at the origin we find 
t — 2/(2 + i[cot(fcL) 4- A/fc]), and therefore the transmis- 
sion through the main wire is 



T= |tp = 4/(4+[cot(/sL) + A/fc]2). 



(21) 



We note that this expression does not reproduce T = 1 
for L = and A = 0, probably due to problems with the 
order of limits. For L > 0, T oscillates with k between 
(full reflection) and 1 (full transmission). Sample results 
are shown in Fig. [21 Since T(fc) = T{—k), we show only 
results for A; > 0. When A = 0, T(fc) is periodic, with 
period tt/L. However, for a non-zero energy-independent 
A the transmission is no longer periodic, approaching the 
zero-A behavior only for large k. 




FIG. 2; Transmission for the continuum stub, 
(dashed) line represents A = (A = 3). 



The full 



We next consider the same model using the discrete 
equations. With a finite lattice constant a, the wave 
function is written as in Eq. (|19p . with x replaced by na. 
Similarly, we set L = {N + l)a and replace y in Eq. pO)) 
by ma, with m — 0,1, ..., A^ -|- 1, where N represents the 
number of bonds on the stub (5 in Fig. [T]), after setting 
tjJs{N -I- 1) = 0. The tight binding equation at the origin 
then yields 



t = 2/ {2 + i[cot[ka{N + 1)] + eo/[Jsin(A:a)] 



cot(A:a)]), 
(22) 



and T = |<p. It is easy to sec the replacement given in 
Eq. HH). Note that (unUke Eq. ^, Eq. ^ does 
reproduce T = 1 for TV = and eo = 0. 

Figure [3] compares the continuum resuh with A = 
(same as in Fig. [5]) with the tight binding result with cq = 
0, for a = L/5. As seen from Eq. (|22p . the transmission 
is now periodic with the new larger period n/a (= Stt/L 
for the case drawn in Fig. [3]). Therefore, we show only 
the lower half of the energy band, < A: < ir/a. As 
expected, the two models differ mostly near the edges 
of the energy band, where cot(fca) diverges. In contrast, 
the two approaches give very close results near the center 
of the band, e ^ or ka ^ tt/2. However, the difference 
between the two models increases as one moves away from 
this region. 




FIG. 3: Transmission for the continuum stub (full line), with 
A = 0, and the tight binding stub (dashed line), with eo — 0, 
ioi N + l = L/a = 5. 

The above results obviously change with the parame- 
ters A and eo. For example. Fig. 3] compares the contin- 
uum model with A = with the tight binding model with 
^o/J = (A/o — 2) = 1, for several values of a. As might 
be anticipated from Eq. p4|) , this value of eg brings the 
tight binding equations close to the continuum ones near 
the bottom of the band, i. e. at small k. Indeed, Fig. 
m shows almost a coincidence of the results near fc = 0. 
This coincidence improves as TV -I- 1 ~ L/a increases. 
However, the results are very different near the center of 
the band, ka ^ 7r/2. 



IV. Conclusions 

A detailed comparison between the tight binding 
"atomic" approach and the continuum approach to quan- 
tum networks reveals that the two models are generally 
different from each other, and one cannot expect them 
to yield the same results for all energies. However, the 
two models agree quantitatively near the center of the ID 
band (which describes each of the wires in the network) , 
e '^ 0. The two models are definitely very different near 
the band edges, where one encounters the van Hove sin- 
gularities. Some of the differences between the models 
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FIG. 4: Same as Fig. [3 but with eo/J = Mj - 2 = 1. The 
dashed curves correspond to A'^-f 1 = 5, 10 and 20 (increasing 
dashes) . 



can be removed by adding appropriate site energies to 
the tight binding model, or by tuning the parameters of 
the matching conditions at the vertices of the continuous 
graphs. Sending the lattice constant to zero brings the 
two models close to each other near the bottom of the 
band only for very special choices of the site energies. 

This paper presented only one simple example of a 
scattering solution, to demonstrate our main points. In 
such a scattering solution, where one is free to choose 
the energy of the scattered electrons near e = 0, the dif- 
ference between the two approaches is small. This free- 
dom does not exist in other situations. Foe example, 
we have recently followed Bcrcioux et alM-'^ in consid- 
ering an infinite chain of diamonds, where the electrons 
are subject to both the Aharonov-Bohm flu5*2^ and the 
Rashba spin-orbit interactions!^ With these two effects, 
the various wave functions become 2-component spinors, 
and the basic hopping terms become 2x2 unitary ma- 
trices, which also contain phase factors. Furthermore, 
the probability current contains additional terms, due to 
these phases. Unlike Bercioux et al., who found the spec- 
trum of the diamond chain within the continuum model, 
we solved the same model within the tight binding ap- 
proach, where each edge of each diamond represented 
a single bond.— We have since repeated this calcula- 
tion with more "atoms" within each edge. Interestingly, 
the qualitative features of the resulting energy bands in 
the two approaches are similar. However, the interest- 
ing physics (e. g. the spin filtering found in Ref. 1271 ) 
does involve moving away from the center of the basic 
ID band. Therefore, the two approaches do yield differ- 
ent quantitative results. The continuum and the discrete 
approaches certainly give different results when one cal- 
culates the spectrum of a finite network. In this case, the 
resulting energies can assume any value, including values 
which are outside the "basic" band. In the latter case, 
the related fc's are imaginary. 

In the above discussion we have not attempted to judge 
which approach is "better" . As we mentioned in pass- 
ing, both approaches are sometimes oversimplified, miss- 



ing some important physics. In practical situations, it is 
probably best to use both approaches and compare the 
results. 
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